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STABILITY OF A FLOWING SYSTEM USING INFORMATION THEORY 

by N. Stankiewicz 
Lewis Research Center 

SUMMARY 

The principle of maximum entropy (disorder) is extended to the nonequilibrium 
situation of viscous flow and a condition for ’’thermodynamic” stability is found. Vio- 
lation of this condition coincides with boundary layer separation in a Newtonian fluid. 


INTRODUCTION 

The concept of maximum entropy is well known in classical thermodynamics, where 
it is used as the definition of the stable equilibrium condition. However, information 
theory with its assumption that the disorder of a statistical system is a measure of its 
entropy can provide a way of extending the maximum entropy (maximum disorder) con- 
cept to nonequilibrium situtations. 

The application of information theory to a system of particles having only continuous 
translational degrees of freedom leads to a thermodynamic description of an ideal, 
classical gas. 

For a compound gas including internal degrees of freedom, the method lends itself 
easily, for example, to the nonequilibrium phenomenon of the two temperature behavior 
of a plasma (electron heating). In reference 1 this approach is applied to a plasma 
where the electrons, ions, and neutrals (with internal energy levels) are treated as 
reacting species. For the case of a quiescent plasma, the Saha equation at the electron 
temperature holds as expected; however, for the case of a current-carrying plasma, a 
departure from ’’Saha equilibrium” occurs. This effect is negligible at ordinary 
currents but might be important in high -density areas or in some shock phenomena. 

The variational method used in information theory directly yields the most probable 
distribution function consistent with the constraints imposed on the system. In applying 
the technique to a flowing system, the constraints are the various flow moments of the 
distribution function. This is in contrast to equilibrium systems where only static con- 
straints are used. The constraints are not numerical constants but may vary from 



point to point as determined from the hierarchy of flow equations. These flow equations 
are derived by taking moments of the Boltzmann collision equation. Unfortunately this 
gives an open set in which each equation involves the next higher flow moment. 

In reference 2 the relation between the most probable velocity distribution function 
and the flow equations is investigated for a specific case in which the system constraints 
are mass density, momentum density, and stress energy density. Because the form of 
the distribution function is known from maximizing the entropy, all the required mo- 
ments can be calculated, thus closing the set of flow equations . For the special case of 
one -dimensional channel flow, the equations give Newton's law of viscosity and a para- 
bolic flow profile. The distribution function (the ellipsoidal velocity distribution func- 
tion) is therefore shown to be consistent with a system undergoing viscous flow. 

The stability of a system is a natural consideration in equilibrium thermodynamics 
and also in the study of flowing systems. In equilibrium thermodynamics there is a set 
of stability conditions (see e. g. , Callen, ref. 3) which is derivable from a minimum 
energy principle. A violation of these conditions indicates that the system is undergoing 
a phase transition. That is, the system separates into two or more subsystems. Such 
a breakup is suggestive of the onset of boundary layer separation in a flowing system. 

This study was undertaken to find the stability conditions for a flowing system using 
the variational principle of maximum entropy and to determine whether a violation of 
these conditions coincides with a physically recognizable event. 

The analysis will be limited to the special case of a system whose particles have an 
ellipsoidal velocity distribution function. This case is consistent with a flowing system 
undergoing viscous dissipation. 


SYMBOLS 

3 

d V dvj dvg dvg 

F phase space density 

f distribution function 

h Planck's constant 

m particle mass 

n particle number density 

S entropy density 

Vj 1 component of particle velocity 

Fjj Lagrange multiplier conjugate to n(VjVj) 
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6 variation of 

Mj Lagrange multiplier conjugate to n<Vj) 

i}/ Legendre transformation of entropy 

0 ) degeneracy 

Lagrange multiplier conji^ate to n 
< > average value of 


ANALYSIS 

The variational equation which gives the ellipsoidal distribution function is 


6S = S25n MjC (n <v.>) + ^ 5 (n (v.v. >) (1) 

i hi 

where Si, and Ty are Lagrange multipliers. The entropy is defined as (ref. 1) 

S = f f - 1) (2) 


and the constraints are 



(3) 


(4) 

n<VjVj> = dv® 

(5) 


The integrals are carried out over all velocities from -«> to +«>. 

The variation 6 in equation (1) operates on the distribution function f . With the 
substitution of equations (2) to (5), this yields the ellipsoidal distribution function 
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( 6 ) 


f = exp 




The relation between f and the density F in phase space is 


F = f (7) 

O 

The maximum density in phase space is w(m/h) for the completely degenerate gas, 
for which case f = 1 . 

Because S is a maximum consistent with the constraints of the system, its second 
order variation must be negative, that is. 


6^S < 0 (8) 

This condition must be satisfied for stability. However, it is convenient mathematically 
to work out the stability conditions using a Legendre transformation of the entropy, 
namely. 


Then 


^ =S - SJn rjjn<VjVj> 

i ij 


(9) 


= -nett n<v.>6Mi 

i iJ 


( 10 ) 


It is a property of Legendre transformations (see ref. 4) to preserve the variational 
principle. Thus the function j// is also an extremal but minimized, that is. 


eV = 0 (11) 

It can be shown, using equations (2) to (6) and the transformation of equation (9) that 

4^ = n (12) 
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Therefore, from equation (10) 



These relations can be verified by carrying out the partial derivatives on equation (3) 
and using the distribution function in equation (6) . 

The second order variation of n (=»//) is then (ref. 3) 





The coefficients of the variations in equation (14) form a 10 by 10 matrix. (There 
are 10 independent Lagrange multipliers: three components of the vector ju, and 

six elements of the symmetric tensor r. ) The quadratic form of equation (14) is 
greater than zero (see eq. (11)), if the determinant of the matrix is positive definite. 
A necessary and sufficient condition for positive definiteness requires that the leading 
principal minors are positive (see e.g. , ref. 5). 

Because the matrix is of order 10, there will be 10 independent conditions for 
stability that can be found. However, the rows and columns of the matrix can be 
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arranged in many ways without disturbii^ the value of its determinant. As -a conse- 
quence, only the principal minors of order one and two need be considered; all other 
minors will give redundant conditions. 

The 10 basic stability conditions are then found to be 




> 0 


( 15 ) 


A 

dlA 


>0; i= 1,2, 3 


(16) 





> 0; i = 1,2,3 


(17) 



3^n 


> 0; Mi = 1,2, 3 


(18) 


Using the definition of n in equation (3) and the distribution function from equation (6) 
and carrying out the partial derivatives, the stability conditions then are 


n>0 


(15a) 


n<vf) >0 


n ^(< Vi > - < Vi >^)>0 


n^(<ViVj> - <Vj) <v.)j 


> 0 


(16a) 

(17a) 

(18a) 


Equations (15a) and (16a) are trivial and are always satisfied. Equation (17a) is a 
condition that the mean random energy for each direction be positive. This is equiv- 
alent to the condition that the temperature be positive and is consistent with the ideal 
gas assumption. For real gases the analogous condition should be that the temperature 
be greater than a temperature at which a change of phase occurs. 

The shear stress is defined in kinetic theory as mn (<Vj^Vj) - (Vj> <Vj>) and must be 
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positive accordii^ to equation (18a) if the flow is to remain stable. For Newtonian 
fluids the shear stress at the wall is proportional to the velocity gradient and a violation 
of equation (18a) will occur whenever the gradient vanishes or changes sign. The 
vanishing of the velocity gradient at the wall is recognized in fluid dynamics as the onset 
of boundary layer separation. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 8, 1970, 

129-02. 
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